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The dependence of transverse momentum spectra of neutral pions and n mesons with pt < 
16 GeV/c and px < 12 GeV/c, respectively, on the centrality of the collision has been measured 
at mid-rapidity by the PHENIX experiment at RHIC in d+Au collisions at ^/s N n = 200 GeV. 
The measured yields are compared to those in p + p collisions at the same ^/s NN scaled by the 
number of underlying nucleon-nucleon collisions in d+Au. At all centralities the yield ratios show 
no suppression, in contrast to the strong suppression seen for central Au+Au collisions at RHIC. 
Only a weak px and centrality dependence can be observed. 
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High-energy nucleus- nucleus collisions provide the op- 
portunity to study strongly interacting matter at very 
high energy densities where Quantum Chromodynamics 
(QCD) predicts a transition from normal nuclear matter 
to a deconfined system of quarks and gluons, the Quark- 
Gluon Plasma (QGP) yj. At the Relativistic Heavy Ion 
Collider (RHIC) the energy density is well in excess of 
the critical energy density that is expected for this tran- 
sition One of the most intriguing results observed 
at RHIC so far is the suppression of hadrons with high 
transverse momentum (px ) m central (head-on) Au+Au 
collisions. The hadron yield at high px is a factor of 5 less 
than expected from p + p collisions scaled by the num- 
ber of corresponding nucleon-nucleon collisions U] . Such 
suppression was predicted as an effect of parton energy 
loss in the medium generated in the collisions 0, 0] . A 
control experiment of d+Au collisions, where no medium 
is produced in the final state of the collision, showed no 
indication of hadron suppression at mid-rapidity [6j , rul- 
ing out strong initial-state effects (final-state energy loss 
in the cold nucleus is generally expected to be small) as 
the cause for the suppression in Au+Au. For a better 
understanding of the medium effects at work in Au+Au, 
however, it is crucial to explore the exact role initial- 
state effects play in the modification of high-px particle 
production at RHIC. 

Initial-state nuclear effects include the Cronin effect, 
shadowing, and gluon saturation. The Cronin effect Q, 
an enhancement of the particle yield at intermediate px > 
is usually attributed to multiple soft parton scatterings 
before a hard interaction of the parton (pj broadening) . 
The shadowing of the structure function [8j modifies the 
particle yield depending on the parton momentum frac- 
tion, XBj, probed in the partonic scattering. An alterna- 
tive model of the initial state of a nucleus is the gluon 
saturation or color glass condensate (CGC) in which 
the gluon population at low XBj is limited by non-linear 
gluon-gluon dynamics. In this picture, particle produc- 
tion at moderate j?x originally was predicted to be sup- 
pressed in central d+Au collisions at RHIC Q. In re- 
cent CGC models, a Cronin enhancement can also be 
reproduced with a suitable choice of initial-state param- 
eters 

One established way to test the contribution of differ- 
ent initial- and final-state nuclear effects is the study of 
the centrality dependence of particle production at high 
px • Initial state and medium effects are strongest in 
central collisions. In Au+Au collisions a strong depen- 
dence of the suppression of high-px hadrons on the cen- 
trality of the collision has been observed 

[HE1E3; the 

suppression weakens going to peripheral collisions and 
finally disappears. This can be compared to the cen- 
trality dependence of (initial-state) hadron production 
in d+Au. The yield of non-identified charged hadrons 
in d+Au collisions with px < 6 GeV/c was found to 
be increasingly enhanced going from peripheral to cen- 



tral collisions 14, mainly attributed to the influence of 
(anti) protons 15j. At high px the baryon contribution 
to the yield of unidentified charged hadrons is expected 
to become small and instead the yield is dominated by 
charged pions All this sparks paramount interest in 
the centrality dependence of neutral pion (jr ) production 
especially as it can be measured up to very high px where 
particle production is truly perturbative. Furthermore, 
the high-px measurement of an additional identified par- 
ticle like the eta meson (77), with four times the mass of 
the pion, may shed light on the question to what extent 
the particle-species dependence of the suppression (en- 
hancement) observed in Au+Au (d+Au) depends on the 
number of constituent quarks rather than on the mass of 
the particle [3, E3] . 

In this Letter we present measurements by the 
PHENIX experiment |l8| on the production of 7r° and 
7] in p + p and d+Au collisions at s NN — 200 GeV. The 
data provides the first measurement of neutral mesons in 
d+Au collisions at mid-rapidity as a function of the cen- 
trality of the collision. The tt° measurements described 
in this paper are similar to the analysis of minimum bias 
d+Au data in 6] but are based on an improved data 
set that allows the study of the particle production for 
different selections of the centrality of the collision. 

7T° and r] are measured by the PHENIX electromag- 
netic calorimeter (EMCal) via the ir° — > 77 and r\ — > 77 
decay. The EMCal consists of six lead scintillator (PbSc) 
and two lead glass (PbGl) sectors, each located at a ra- 
dial distance of ~ 5 m from the beam axis. The detector 
covers a pseudorapidity range of \r\\ < 0.35 and an az- 
imuthal angle of A<fi — it. The EMCal granularity is 
Ar) xA(4k 0.011 x 0.011 for the PbSc and 0.008 x 0.008 
for the PbGl. The data sets from PbSc and PbGl are an- 
alyzed separately and combined for the final results. The 
energy calibration for the EMCal is obtained from beam 
tests, cosmic rays, and minimum ionizing energy peaks 
of charged hadrons. In a recent improvement of the cal- 
ibration, the EMCal is calibrated by the invariant mass 
distribution of neutral pions for each of the 24768 read- 
out channels separately. The uncertainty on the energy 
scale is 1.2%. 

The data used in this analysis was recorded in 2002- 
2003 (RHIC Run-3) under two different trigger condi- 
tions. 25.2 xlO 6 and 58.3 xlO 6 minimum bias events 
were analyzed for p+p and d+Au collisions, respectively. 
Minimum bias (MB) events are triggered by the Beam- 
Beam Counters (BBC) [3 (\r)\ = 3.0-3.9) and require 
a vertex position along the beam axis within \z\ < 30 
cm. The minimum bias trigger accepts (88 ±4)% of all 
inelastic d+Au collisions that satisfy the vertex condi- 
tion. This corresponds to 1.99 b ±5.2%, the measured 
fraction of the total d+Au inelastic cross section, deter- 
mined using photo-dissociation of the deuteron |lj| . In 
p + p this trigger measures 23.0 mb ±9.7% of the p + p 
inelastic cross section. The measured particle yields are 
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corrected for the p + p MB trigger bias 20] : the MB 
trigger measures only (79 ± 2)% of high-px particles. In 
d+Au collisions this fraction varies from 85% to 100% 
from peripheral to central collisions; here the uncertainty 
is ~ 3%. The second data sample was collected with a 
high-pT photon trigger in the EMCal in addition to the 
MB trigger requirement in order to extend the measure- 
ment to higher px ■ This trigger requires a photon of 
p T > 1.4(1.4) and p T > 2.5(3.5) for PbSc (PbGl) and 
for p + p and d+Au collisions, respectively. We analyzed 
45.1 xlO 6 (19.5 xlO 6 ) events in p + p (d+Au) under 
this trigger condition. The sampled integrated luminos- 
ity was 216 nb" 1 for p + p and 1.5 nb _1 for d+Au. (In 
d+Au that corresponds to an integrated nucleon-nucleon 
luminosity of 590 nb _1 ). 

The division of d+Au collisions in different centrality 
classes is based on the charge deposited in the backward 
BBC (—3.9 < rj < —3.0), i. c. in the Au beam direction. 
For each centrality class the corresponding average nu- 
clear overlap function (Tab) is calculated using a Glauber 
Monte Carlo model and simulations of the BBC, taking 
into account its limited efficiency for peripheral collisions. 
For the four centrality classes (0-20%, 20-40%, 40-60% 
and 60 — 88%) used in this analysis, (Tab) translates 
into an average number of nucleon-nucleon collisions per 
A+B collision, (JVccji) = of*, x (Tab) , of (15.4+1.0), 
(10.6 ± 0.7), (7.0 ± 0.6) and (3.1 ± 0.3), respectively. 

Photon candidates in the EMCal are selected by ap- 
plying particle identification (PID) cuts based on the 
shower profile in the detector. To determine the yields 
of 7T° and r], the invariant mass of all photon pairs with 
an energy asymmetry \E\ — E 2 \/{Ei + E 2 ) < 0.7 in a 
given px bin is calculated. After subtraction of the com- 
binatorial background the invariant mass distribution is 
integrated around the particle mass peak [ll|; the in- 
tegration window reflects thereby the px dependence of 
the mass peak position and width. The combinatorial 
background is determined by pairing photons from dif- 
ferent events with similar centrality (for d+Au) and ver- 
tex. In this analysis the signal-to-background ratio for 
high-px tt° is about 25 and 13 at px =4 GeV/c in p + p 
and central d+Au collisions, respectively. It decreases 
to 7 and 2 at px =2 GeV/c. For 77, this ratio is about 
2 at pt —8 GeV/c, decreasing to 0.3 (p + p) and 0.2 
(central d+Au) at px =3 GeV/c. The raw spectra are 
corrected for trigger efficiency, acceptance, and recon- 
struction efficiency. This includes dead areas, the influ- 
ence of energy resolution, analysis cuts, the peak extrac- 
tion window, and photon conversion. The corrections are 
determined using Monte Carlo simulations. Due to the 
fine granularity of the calorimeter, occupancy effects are 
negligible. Furthermore, the ir° spectra are corrected at 
Pt >10 GeV/c (15 GeV/c) for two-photon merging ef- 
fects in the PbSc (PbGl), studied in Monte Carlo sim- 
ulations and confirmed with test beam data. Finally 
a correction in the 7r° and 77 yields to account for the 



true mean value of each px bin is applied to the steeply 
falling spectra. For px < 3.5 (3.0) GeV/c the p + p 7r° 
(77) spectrum is calculated from the minimum bias data 
sample, above this threshold the high-px triggered sam- 
ple is used. In d+Au, this transition is made at px = 4.5 
and 3.5 GeV/c for tt° and 77, respectively. 

The main contributions to the systematic uncertainty 
on the p+p and d+Au spectra are given in Table [I] for 
7T° and 77. Most uncertainties are identical for p + p and 
d+Au, only the uncertainty on the peak extraction is 
slightly larger in d+Au. Category (d) includes uncertain- 
ties on the EMCal global energy scale and non-linearity. 
The uncertainties in (d) and (e) are partially correlated. 
All others are uncorrelated. 

The fully corrected px distributions of ir° and 77 are 
shown in Fig. ^ The top panels show the invariant yield 
in d+Au collisions for four centrality bins scaled for clar- 
ity by the factors indicated. The bottom panels show 
the invariant cross section in p + p and d+Au collisions. 
The improved dataset allows the study of n° (77) produc- 
tion up to 18 (12) GeV/c, the highest px values measured 
for identified particles in p(d)+A collisions. For the first 
time, the invariant cross section for 7r° and 77 in d+Au 
collisions has been measured at this energy. The ir° result 
in p + p agrees with the previous measurement at y/s NN 
= 200 GeV [2fj| within statistical uncertainties, and con- 
firms the agreement with pQCD within the uncertainty 
of the calculation. Therefore the p + p cross section can 
be used as a well-understood reference for the production 
in d+A and Au+Au collisions. 

To quantify nuclear medium effects at high px it is 
customary to use the nuclear modification factor which 
is given by the ratio of the d+Au yield to the p+p cross 
section [ll| scaled by (Tab)'- 



Rab(pt) 



d 2 Nl B /dyd PT 
(Tab) x d 2 a;°/dydp 



(1) 



The average nuclear overlap function (Tab)-, averaged 



TABLE I: Main systematic uncertainties in % on ir° and 77 
spectra. The uncertainties are given for PbSc (PbGl). The 
normalization uncertainties of 9.7% for the p + p and 5.2% 
for the d+Au cross section as well as the MB-trigger-bias 
uncertainty of ~ 3% for the centrality-selected yields are not 
listed. 
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FIG. 1: Top: invariant yields at mid-rapidity for n (left) and 
77 (right) in d+Au collisions as a function of pt for different 
selections of the centrality of the collision. Bottom: invariant 
cross section at mid-rapidity for 7r° (left) and r? (right) in p+p 
and d+Au collisions as a function of pt ■ 
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FIG. 2: Nuclear modification factor R^a for it and 77 in 
different centrality selections and min. bias data. The bands 
around the data points show systematic errors which can vary 
with pt ■ The shaded bands around unity indicate the (Tab) 
uncertainty and the small bands on the left side of the data 
points indicate the normalization uncertainty due to the p + p 
reference. 



over the respective impact parameter range, is deter- 
mined solely by the density distribution of the nuclcons 
in the nuclei A and B and the impact parameter. 

Figure [21 shows the nuclear modification fac- 
tor R<ia{pt) for 7T° and 77 in d+Au collisions at 
x/s NN —200 GeV for four different centrality selections 
and for minimum bias events. As the p + p and d+Au 
measurements were both made in the same year, many 
of the systematic errors associated with detector perfor- 
mance were nearly identical and the corresponding sys- 
tematic errors in the comparison are negligible. Within 
systematic errors RdA(PT) for 7r° and 77 is « 1 in all 
centrality bins, and only a weak px dependence can be 
seen. In order to check the absolute normalization sys- 
tematics we can also calculate RdA^Pr) using the inelas- 
tic cross section measured through photo-dissociation of 
the deuteron. This constitutes an important cross check. 
It replaces the systematic uncertainties of the BBC effi- 
ciency and (Tab), which are determined by model cal- 
culations, by the uncertainty of the cross section mea- 
surement of similar size. The resulting RcLa{pt) is 9.8 % 
larger than that obtained from the minimum bias yield, 
consistent within 1.5 a. 



Though very different in mass, 77 and 7r° show a simi- 
lar, weak centrality dependence of RdAijpT) over the mea- 
sured pt range. These results do not show the significant 
enhancement seen for protons where the proton Raa is 
substantially larger than that of pions in the intermedi- 
ate p T (2 GeV/c < p T < 4 GeV/c) region 0. The tt° 
data exhibit small shape variations with centrality that 
may be due to initial-state effects including shadowing 
and multiple scattering. Possible Cronin enhancements 
in the intermediate pt region due to initial-state multi- 
ple scattering or anti-shadowing are not more than 10% 
around 4 GeV/c. At low pt (pt < 3 GeV/c) the drop to- 
wards smaller RdA is consistent with analogous measure- 
ments for charged pions [15j and is usually attributed 
to a change to a regime of soft physics (N part scaling) 
at the smallest pt values. At the largest pt values mea- 
sured (px > 9 GeV/c) the most central 7r° result hints at 
a small suppression, though this is only a ~ 1.7 sigma 
effect. 

In conclusion, we have presented the first study of the 
centrality dependence of 7T° and r\ production at mid- 
rapidity in g?+Au collisions at s NN = 200 GeV. Trans- 
verse momentum spectra up to pt = 18 and 12 GeV /c 
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have been measured for 7r° and 77, respectively. The in- 
variant yield per nucleon-nucleon collision is compared to 
that in p + p collisions measured at the same y/s NN . The 
strong suppression observed for 7r° production at high 
Pt in central Au-Au collisions cannot be seen for d+Au in 
any centrality: Within systematic errors RcIa(pt) is » 1 
in all centrality bins. A weak centrality dependence of the 
shape of RdA versus pi is seen, presumably due to initial- 
state effects. A possible Cronin enhancement is substan- 
tially smaller than the i?<M>1.9 that corresonds to re- 
sults from lower energy measurements 0, El- Within 
systematic errors RdA for and r\ agrees well, giving 
no indication for cold nuclear matter effects having a 
mass dependence. Since nuclear modifications in d+Au 
are small even in the most central collisions where initial 
state effects are expected to be largest, we conclude that 
initial state effects in Au+Au must be small as well, so 
the large suppression seen in Au+Au must be mostly due 
to medium effects. 

We thank the staff of the Collider- Accelerator and 
Physics Departments at BNL for their vital contri- 
butions. We acknowledge support from the Depart- 
ment of Energy and NSF (U.S.A.), MEXT and JSPS 
(Japan), CNPq and FAPESP (Brazil), NSFC (China), 
IN2P3/CNRS, CEA, and ARMINES (France), BMBF, 
DAAD, and AvH (Germany), OTKA (Hungary), DAE 
and DST (India), ISF (Israel), KRF and KOSEF (Ko- 
rea), RMIST, RAS, and RMAE (Russia), VR and KAW 
(Sweden), U.S. CRDF for the FSU, US-Hungarian NSF- 
OTKA-MTA, and US-Israel BSF. 



* PHENIX Spokesperson: zajc@nevis.columbia.edu 



' Deceased 

[1] J. W. Harris and B. Miiller, Ann. Rev. Nucl. Part. Sci. 
46, 71 (1996). 

[2] K. Adcox et al. (PHENIX), Nucl. Phys. A757, 184 
(2005). 

[3] K. Adcox et al. (PHENIX), Phys. Rev. Lett. 88, 022301 

(2002) . 

[4] J. D. Bjorken (1982), FERMILAB-PUB-82-059-THY. 
[5] R. Baier, D. Schiff, and B. G. Zakharov, Ann. Rev. Nucl. 

Part. Sci. 50, 37 (2000). 
[6] S. S. Adler et al. (PHENIX), Phys. Rev. Lett. 91, 072303 

(2003) . 

[7] J. W. Cronin et al., Phys. Rev. Dll, 3105 (1975). 
[8] M. Arneodo, Phys. Rept. 240, 301 (1994). 
[9] D. Kharzeev, E. Levin, and L. McLerran, Phys. Lett. 
B561, 93 (2003). 
[10] D. Kharzeev, Y. V. Kovchegov, and K. Tuchin, Phys. 

Rev. D68, 094013 (2003). 
[11] S. S. Adler et al. (PHENIX), Phys. Rev. Lett. 91, 072301 

(2003) . 

[12] S. S. Adler et al. (PHENIX), Phys. Rev. C69, 034910 

(2004) . 

[13] J. Adams et al. (STAR), Phys. Rev. Lett. 91, 172302 

(2003) . 

[14] S. S. Adler et al. (PHENIX), High p T charged central- 
ity dependence and h/ir° ratio in d+Au collisions, to be 
published. 

[15] S. S. Adler et al. (PHENIX) (2006), nucl-ex/0603010. 
[16] S. S. Adler et al. (PHENIX), Phys. Rev. C72, 014903 

(2005) . 

[17] R. C. Hwa and C. B. Yang, Phys. Rev. Lett. 93, 082302 

(2004) . 

[18] K. Adcox et al. (PHENIX), Nucl. Instmm. Meth. A499, 
469 (2003). 

[19] S. N. White, AIP Conf. Proc. 792, 527 (2005). 
[20] S. S. Adler et al. (PHENIX), Phys. Rev. Lett. 91, 241803 
(2003). 

[21] A. L. S. Angelis et al. (BCMOR), Phys. Lett. B185, 213 
(1987). 



